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1. Introduction 

There is by now strong and growing evidence that the very early universe went through a 
phase of inflation [1]. What led to this phase of accelerating expansion, and what the un- 
derlying microscopic theory is, are questions that are still under intense investigation. Our 
current understanding from the nearly scale invariant spectrum of the cosmic microwave 
background (CMB) anisotropics is that such an inflationary phase (at least the observable 
e-folds) can be well described by an approximate de Sitter universe. 

An interesting modification of a (quasi-) de Sitter spacetime that would still be con- 
sistent with an inflationary history of the universe is the introduction of some primordial 
radiation. At late times, this radiation would get inflated away and the naive expectation is 
that this radiation will therefore not leave any observable effects on the CMB temperature 
anisotropy power spectrum. Various cosmological models have been proposed in which it 
is natural to include some amount of radiation in the early (quasi-) de Sitter phase [2-4] 
and some general results on radiation-filled FRW universes can be found in [5]. Moreover, 
going beyond the mini-superspace approximation in selecting the wave function of the uni- 
verse [6-8] naturally leads to an effective description in terms of a radiation deformed de 
Sitter geometry, where the radiation density is related to the fundamental (string) length 
scale in the problem [9,10]. 1 

1 A similar conclusion was reached in [11] and [12] based on string thermodynamic arguments and quan- 
tum gravitational loop effects, respectively. 
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All these considerations warrant a more detailed study of a radiation deformed de 
Sitter geometry. In particular, it would be interesting to see whether primordial radiation 
could leave some observational signatures. Here we report that in a de Sitter plus radiation 
universe, remarkably, the power spectrum of scalar field fluctuations is determined by a 
known differential equation: Lame's equation. This allows for a precise analysis of the 
spectrum and the modifications with respect to a pure de Sitter geometry. 

Scalar fields in pure and deformed de Sitter spacetime have of course been well studied. 
[14,15] studied the mode functions of a scalar field on a pure de Sitter background. The open 
radiation deformed de Sitter geometry is considered in [17]. We study a closed radiation 
deformed de Sitter geometry. There is a crucial difference between the open and the closed 
geometries. In an open radiation deformed geometry, the universe begins in a singularity, 
followed by a phase of radiation domination, after which the universe becomes dominated 
by the cosmological constant. In contrast, the radiation deformed global de Sitter spacetime 
interpolates between a pure closed de Sitter and an Einstein static universe as the total 
amount of radiation is varied from zero to the maximally allowed value, and there is no 
singularity in the geometry. The presence of any more radiation (than the maximum 
allowed amount) leads to a spacelike singularity. In Euclidean signature, the geometry 
corresponds to a deformation of an S A to a "barrel" shaped geometry [9,10]. So below the 
maximally allowed value for the radiation density the initial singularity is absent and the 
geometry can safely be analytically continued to a bouncing Lorentzian geometry. 

Our study is motivated in part because the Euclidean geometry allows an interpretation 
as the instanton process involved in selecting the wave function of the universe [6-8]. 
What might be a natural choice for the wave function of the universe is a subject of 
long debate [13], which we will not address in this paper (see [9,10] for discussion). In the 
original Hartle-Hawking context, it allows for an unambiguous determination of the vacuum 
state of fluctuations [23,24]. The cosmological vacuum ambiguity has recently received 
renewed interest as future CMB experiments are sensitive enough to probe it [30]. We will 
point out that in a Hartle-Hawking cosmological scenario with radiation the Bunch-Davies 
initial state is no longer unambiguously selected by analytic continuation to the Euclidean 
geometry. 

To briefly summarize our main result: the effect of the radiation is to lead to a moder- 
ation of the closed dS red tilt of the scalar field power spectrum. The obvious explanation 
for this is that the low /c-modes are sensitive to the spatial curvature, i.e. the compactness, 
of the closed de Sitter geometry. The radiation deformation effectively shrinks the volume 
of the spatial 5 3 's, which changes the red tilt in the fluctuations. Consistency with the 
observed flatness of our universe, however, make the red-tilt due to curvature barely ob- 
servable. Similarly, the additional red-tilt due a radiation deformation is unlikely to show 
up in future experiments. 

2. A closed universe with a positive cosmological constant and radiation 

The scale factor for a closed universe with a positive cosmological constant A plus a radi- 
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ation term with energy density p*/a obeys the equation 

1 d 2 a(t) 1 (da(t)\ 2 1 g. 
a(t) dt 2 + a(t) 2 V dt J a(ty a(t) 4 ' { ' 

We use i for time in Lorentzian signature metric, r for time in Euclidean signature metric, 
natural units in which 8irG = c = H = 1. The solution to Eqn. (|2.l| ) is given by 



1 



a(i) = ^{l + y/l-p cosh(2Lt) J . (2.2) 



Here /3 = 4L 2 /?* and L = y is the late-time horizon (see appendix [B]). The Euclidean 
version is a deformation of S* 4 . We will refer to it as the "barrel" geometry [10] 



a(T) = V2L V V + vWcos(2Lt)J • ( 2 - 3 ) 

As /3 — > (no radiation) one recovers the pure de Sitter geometry. In the opposite limit 
(5 — ► 1 one obtains the Einstein Static universe. The Euclidean geometry corresponds to 
an S A for f3 = 0, a cylinder (i? 1 x S 3 ) for 0=1, and for general /? interpolates between 
these two, resembling a barrel (Fig.(l)). 

Perhaps a geometrically more intuitive way to write the scale factor is 

1 



o(r) = a/A 2 +cos 2 (Lr) . (2.4) 

V ^ L\/1 + 2A 2V V ' ' 

where A 2 = ^^3. In the limit of vanishing radiation, A — > 0, one manifestly recovers 
pure de Sitter. 

Compared to Lorentzian de Sitter (Eq.(p.2j)), the size of the spatial S 3, s shrink as a 
consequence of the radiation, which is most pronounced around the waist at t = 0. At this 
bounce, the waist is reduced by a factor yj ■ This effect reaches a minimum in the 
(critical) limit /3 — >• 1, or equivalently A 2 — > oo. The geometry reduces to a pure de Sitter 
spacetime as t — > ±oo due to the redshifting of the radiation. The Penrose diagram of 
the radiation deformed solution corresponds to a slightly elongated version of the perfect 
square of pure de Sitter [18]. This shows that the compact nature of the spatial slices can 
in principle be ascertained in a finite amount of proper time by fiducial observers, which 
is impossible in pure de Sitter. 

3. Scalar field in a closed de Sitter universe with radiation 

To show that the fluctuation spectrum is determined by the Lame equation, consider a 
minimally coupled scalar field 

1 M 2 

U = --g ab (d a <S>)(d b <S>)-—$ 2 . (3.1) 
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Figure 1: Radiation changes the Euclidean S 4 to a Euclidean Barrel. 
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with radiation 
(P*0) 

> a(t) 



pure de Sitter 
(P = 0) 



Figure 2: The solid figure corresponds to the pure de Sitter universe (Euclidean S 4 ). The dotted 
figure corresponds to a radiation deformed de Sitter universe (Euclidean Barrel). More radiation 
leads to more shrinking at the waist. 



The wave equation (in Euclidean signature) is given by 

A 3 



:i <MlL 

dr 2 a(r) dr a(r) 2 



+ M z U> = , 



(3.2) 



where a = da(r)/dr, and A3 is the Laplacian on the unit three-sphere. As is standard, we 
decompose the field into S 3 spherical harmonics 



$ = fk(r)Y Mm (x,0, t 



(3.3) 



with eigenvalues 



A 3 y fc , m = -(A;-l)(A ; + l)y, 



klm 



(k- 1) > / > \m\ > 0. 



(3.4) 
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The mode function must then obey the following equation 



fk — 3-A + 
a 



(k 2 ~ 1) , r2 



a 2 



fk = 0. (3.5) 



We are following the common convention in a closed cosmology to let the minimum 
(co- moving) momentum variable k equal 1, corresponding to the homogeneous mode of 
Eq.(|3.4). Note that in our notation the scale factor a(r) has dimensions of length, imply- 
ing that the co-moving momentum k = a(r)p is a dimensionless number. An important 
thing to realize is that in these conventions the same co-moving scale in initially different 
closed universes (for example one with or without primordial radiation) , will not correspond 
to the same physical scale at late cosmological time. 

In what follows, we shall simplify the above equation and show that the mode functions 
//t(r) are given in terms of solutions to Lame's equation. First make the substitution 

F k (r) = o(t)A(t) . (3.6) 

then Fk must satisfy 

Fk + ^F k -(* + %)Fk-( ( ^± + M 2 )F k = V. (3.7) 
a \a a z J \ a z J 

From the trace of the Einstein equations it is easy to show that for a closed FRW universe 
(in Euclidean signature) 

a a 2 \ 9 1 , 

- + -o ) = "2£ 2 + -y • 3.8 



This leads to the equation 

F k + ^F k + L l Ufa + 1) - ^1 F k = , (3.9) 

where /i(/x + 1) = 2 — M 2 /L 2 . Substituting the explicit form of the scale factor, we obtain 
the following equation 

! + *) F k + Lmt( Lr Wt + L 2 U + 1) + A2( ^ +1)) 2 - t ; (1 + 2A2) ) F k _ 
sin {Lt ) J V sm (Lr) / 

(3.10) 

in terms of f = r + ir/2L. In the limit A — > 0, one recovers the associated Legendre 
equation for pure de Sitter space using global coordinates (closed S 3 slicing) [19,20]. 

Let us now express everything in terms of a variable y = cos 2 (Lf). This yields 



(//(/x + 1) - A: 2 ) + A 2 ( M (^ + 1) - 2fc 2 ) - MQi + l)y 
+ % (y -l)( y -,) Fk{y) = ° ' (3 -^ 
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where z = 1 + A 2 . Eg. (gig) is an ordinary differential equation with four singular points 



(at 0, 1, z,oo). This equation is well studied and is known as Heun's equation [21]. In 
fact, Eq.( 3.11[ ) is a special type of Heun's equation called Lames equation. That scalar 



field fluctuations in closed de Sitter universes plus radiation reduce to Lame's equation 
was also noted in [5]. Here we use this to compute an exact power-spectrum of scalar 
field fluctuations. In this context it is interesting to note that according to [21], Lame's 
equation is the most general equation of Heun type which can arise from a Laplacian after 
separation of variables. For the sake of completeness the general features of Heun's and 
Lame's equation are discussed in appendix [A|. Solutions to Lame's equation exist both as 
a power series and as a series in hypergeometric functions. 

3.1 Initial conditions and vacuum states 

To determine the complete solution to the mode equation ( |3.5|) we must supply initial 
conditions or equivalently define the vacuum state. This is where physics enters. To 
construct the vacuum states we follow [15,16]. In the previous sections we wrote the scalar 
field equation of motion in a metric using Euclidean time. In this section we will work in 



Lorentzian signature, a(t) = Ly /^ +2 ^i y ^ 2 + cosh 2 (Li). 

In order to follow the notation of [16], we start from the Lorentzian counterpart of 
Eq.(p.8|), and define a new time variable r\ given by 

T] = 2arctan (e Lt ) . (3.12) 

As — oo < t < oo, < rj < 7r. For pure de Sitter this variable simply corresponds to 
conformal time, but this is not true in general. In terms of this variable the equation of 
motion becomes 

d 2 F k A 2 sin 77 cos t? ^ / (1 + 2A 2 )fc 2 (2-M 2 /L 2 ) ^ 

+ , , a9 ■ 2 ZT + - , *■>■.- 2 ]*k = U- (3.13) 



dr\ 2 1 + A 2 sin 2 77 drj \ 1 + A 2 sin 2 rj sin 2 77 

It will be useful to recast this equation one more time into a deformation of the defining 
associated Legendre equation by introducing a new function G k related to F k as F k = 
(\/sin 77/(1 + A 2 sin 2 77) 1 / 4 ) G k . Then in terms of the variable x = — cos(rj) one obtains 

(1 _ x2) d^ _ ^ + ((1 + 2A2)(fc2 _ m _ 

v 2 (l + A 2 (l - x 2 )) , 3A 2 (1 + A 2 ) (1 - x 2 ) \ rf . n 

(3.14) 



1-x 2 + 4 (l + A 2 (l-x 2 ))J G ' fc(:E) " ' 



where v 2 = 9/4 - M 2 /L 2 . 

It is worth briefly recapitulating the construction of vacua for de Sitter spacetime 
before we complicate matters by taking into account the radiation term. When A 2 = 0, the 
above equation is by construction an associated Legendre equation, and the time variable 
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r\ corresponds to conformal time. The most general solution is a linear combination of the 
associated Legendre polynomials 

F k ( v ) = y/^j[A k P£_ 1/2 {-coBTi) + B k Q v k _ 1/2 (-cosr))] . (3.15) 

The coefficients A k and B k are not independent. The orthonormality of the physical mode 
functions requires that 

"tHtsH- < 3 - 16 > 

and this leads to the constraint 

A k B* k - B k A% = - r(fc + l/2-^) (3 } 

k k T(k + 1/2 + v) y 1 



A priori there is no preferred vacuum state and all solutions to equation ( 3.17] ) are 



allowed. On physical grounds however, one state, the Bunch-Davies or Euclidean vacuum 
is usually preferred. The Euclidean vacuum is the unique set of initial conditions for which 
the analytically continued Euclidean Green's function has no singularities at the antipodal 
point [22]. For our purposes this can simply be interpreted to mean that the state should 
reduce to the usual flat space vacuum at physical momentum scales much larger than the 
curvature scale, i.e. k/a S> H. This requirement singles out the following values for the 
coefficients [15, 16] 

Ak ~ \AT(k + 1/2 + v) J 

B k = —A k . (3.18) 

In the Hartle-Hawking "tunneling from nothing" scenario for the wavefunction of the uni- 
verse, analyticity has been argued to uniquely resolve the cosmological vacuum ambiguity 
in favor of this Euclidean Bunch-Davies vacuum [23,24]. Clearly this analyticity argu- 
ment no longer holds with even a small amount of radiation present; there is no regularity 
condition at the south pole of the "barrel" . 

3.2 Vacua for the radiation deformed geometry 

We would like to see how this de Sitter result - the one complex parameter family of vacua 
given by Eq.( |3.17| ) - gets modified by the presence of radiation. We do so by a power 
series analysis of solutions to Eq. fl3.14p . Let us denote the two independent solutions 
by G ! fe_ 1 / 2 an d ^&_i/2- We will require that in the limit of vanishing radiation G v k _ 1 i 2 
reduces to P k _i/ 2 (x) and H k _ 1 ^ 2 reduces to Q v k _ X j 2 - Using the power series ansatz G(x) = 
YlrLo c r(^ ~ x) r+a it is straightforward to verify that the two independent solutions to 
Eq.( p,14j ) are given by 

a = v/2, or a = —v/2 , 

£i = 2a(2« + l)-2(fc-l/2)(fc + l/2) ^ - 4a 2 - (k - l/2)(k + 1/2) 

c {2a + 2) 2 -v 2 {2a + 2) 2 -v 2 ' { ' ' 
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Here we have listed only the ratio of the first two coefficients. This suffices to study 
the late time behavior {x — > 1 limit). The coefficient cq is fixed by requiring the so- 
lutions to reduce to the associated Legendre functions for zero radiation. This gives 
Co = ((— l) v T(k + v + 1/2)) / (2 v / 2 v \T(k — v + 1/2)) . The general n-th term recursion re- 
lation is easy to find. 

The above two values of a (viz. a = +v/2, —v/2) give two independent solutions 
an d These are the analogs of P k _ x / 2 an d P^-X/v res P ec ti ver y- To find 

' ^ e analog of Q k _i/% " we remind the reader of the relationship for associated 
Legendre functions [25] 

= TO (^v.w ™(»-> - iro^ -V*)) ■ < 3 - 20 > 

We can find the corresponding relation for H k _ 1 ^ 2 as follows. In Eq.( |3.20|) , the coefficients 
of the associated Legendre functions are independent of the variable x, representing time. 
We expect the same to hold true for the modified relationship. If we can fix these coefficients 
at one instant of time, they will have those values at all times. Now note that at late times 
(as x — ► 1), Eq (|3.14|) takes the form 

(1 - - 2x^M + ((1 + - 1/4) - ^ ) ,(,) = . 

(3.21) 

This is almost the equation for de Sitter spacetime, except for the factor of (1 + 2A 2 ). This 
might seem puzzling. At late times one expects the radiation to have no effect. However, 
recall that even at (fixed) late times, the overall size of a compact spatial slice is scaled down 
due to the presence of the radiation. Now we can construct a relationship analogous to 
Eq.(fj2(J. At late times, G v k _ 1/2 approaches J£_ 1/2 , with ( X 2 - 1/4) = (l + 2A 2 )(/c 2 - 1/4). 
This suggests to simply replace k by x, giving 

= to K'/sM c ° s( ™> - lt-llm G ^ {x) ) ■ (3 ' 22) 

In terms of these two independent solutions G k _^, 2 and i?^_ 1 ^ 2 , the general solution 
to the mode functions reads 

„ ^ V^sinr?) 



(1 + A 2 sin 2 r/) 1 ^ 
Y / smr/ 



G fe (3.23) 
[A k G v k _ 1/2 (- cost]) + B k H^_ 1/2 (- cost])] . 



(1 + A 2 sin 2 7]) 1 / 4 

A careful analysis of the orthonormality constraints shows that the two complex parameters 
A k and B k must be related as 

A k Bl -B k A* k = i y/l + 2A 2 ^(x + 1/2-^) ^ ^ 

fe fc r(x + i/2 + t>) v ; 



The appearance of another \/l + 2 A 2 factor in this condition is important. It appears 
because in the radiation deformed geometry the ry variable does not correspond to conformal 
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time. From Eq^^) one sees that in the late time limit, the scale factor reduces to a(t) ~ 
W1+2A 2 ' ^ S com P are d to pure de Sitter space, this simply amounts to a rescaling of the 
conformal time variable with \/l + 2A 2 . This explains the factor in Eq. ( |3.24 ) . In addition 



we simply replaced k by \ on the right-hand-side of this equation, following our earlier 



argument that at late times the Eq. (|3.14j) reduces to the associated Legendre equation and 
that the normalization condition is independent of time. 

To find the analog of the Bunch-Davies vacuum in the (late time) radiation deformed 
geometry, we proceed in exactly the same way as before. We demand that our solution 
reduces to the usual flat space vacuum at high momentum p = k/a. In addition it should 
reduce to the usual de Sitter Bunch-Davies vacuum in the limit A 2 — ► 0. Keeping in mind 
that (x 2 — 1/4) = (1 + 2A 2 )(k 2 — 1/4), this singles out the following solution 

B k BD = —A k BD . (3.25) 

So we see that the only difference as compared to the pure de Sitter case is the replacement 
of k with x an d the additional factor of (l + 2 A 2 ) 1//4 due to the change in the normalization 
condition Eq.( |3.24D . 

Fundamentally there is deeper difference between the Bunch-Davies in the radiation 
deformed case and that of the pure de Sitter. Were one to insist on analytic behavior 
of the wave function in the Euclidean past, as e.g. in the Hartle-Hawking scenario, then 
one must choose the Euclidean or Bunch Davies vacuum in pure de Sitter [23,24] to avoid 
a singularity at the south pole of the sphere. Because even a tiny amount of radiation 
opens up the South pole to a barrel, there appears to be no reason to insist on well-defined 
behavior at the edge. The Bunch-Davies state is now only preferred for the usual reasons. 

Finally note that one should no longer consider the limit A 2 — ► oo in the late-time 
de Sitter solution; the late-time de Sitter limit is implicitly defined as t > tf, with 
cosh 2 (Ltf) 3> A 2 . Instead we know that the homogeneous solution reduces to the Ein- 
stein Static Universe in the limit A 2 — ► oo. It is not difficult to see that in that case 
the mode-functions will reduce to those of ordinary flat space, i.e. plane waves. As we 
will be interested in a (late-time) inflationary stage, the late-time de Sitter solution is the 
appropriate one to study the leading effects due to the radiation deformation. 



4. The power spectrum of scalar field perturbations 

The large scale density perturbations present in the CMBR temperature anisotropy we see 
today crossed the horizon about 55 e-folds before the end of the inflation. 2 Subsequently 
they remained frozen till horizon re-entry. Assuming we can construct a suitable closed in- 
flationary model based on the radiation deformed de Sitter geometry, the question whether 
there will be any imprint on these perturbations probably depends on whether the radiation 

2 The precise e-folds at which the observable modes crossed the horizon depend on the reheating temper- 
ature. In some extreme cases, the horizon crossing can be as late as 25 e-folds before the end of inflation. 
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was a sizeable component at horizon exit. So we expect that the large wavelength/small 
k perturbations are more likely to be affected by the radiation (as they cross the horizon 
and get frozen in early on during inflation). 

Let us remind the reader that we will calculate the power spectrum of massless scalar 
field perturbations, which is not the same as the primordial power spectrum of scalar 
density perturbations. Nevertheless, as is usual, we will implicitly assume that parts of 
our results continue to hold for the primordial spectrum of scalar density perturbations in 
a slow-roll inflationary generalization of the background geometry, which is the relevant 
quantity responsible for the observable CMBR anisotropy spectrum. 

The advantage of having an exact (late-time) solution to the equations of motion is 
that the power spectrum of fluctuations is precisely given by 

leading term 1,3 

where the correctly normalized mode function for a massless field (v = 3/2) equals 

In the late time limit, just like in the pure de Sitter case, the second term will have 
vanishing contributions, H^^^ — > 0. Using that in the late time limit G^{^ 2 approaches 

the associated Legendre polynomial -P x _ 1/ / 2 ( x )> with asymptotic behavior 

lim IZ(x) = (1 - X )-^-^— , (4.3) 

x~ >1 1.(1 — V) 



the late time mode function is 



IMM^ 2 ) 3/2 ^ + ---- <"> 

To evaluate x, recall that (x 2 - 1/4) = (1 + 2A 2 )(fc 2 - 1/4). Solving for x we get 



V ■ \ 1+2A 2 ^ fc2 V4) . (4.5) 



Thus the exact scalar field power spectrum is given by 



ft_^| A p_(l + 2A>r^^. (4.6, 



This is our main result. For large wavenumber k it reduces to 

+ ...) . (4.7) 



L 2 / 1 



1 + fA 2 
1 + 2A 2 



The overall term corresponds to the well-known de Sitter result. Notice the cancelation 
of the radiation dependent factors between the norm of the mode-function and the explicit 
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X dependence. The term in between brackets corresponds to a small red-tilt (more power at 
small k) of the spectrum of scalar field perturbations. At low values of k, the departure from 
scale invariance is pronounced and the radiation component is of the same order (in 1/k) 
as the effect due to curvature alone [19]. The red tilt is of course easy to understand. Low 
fc-values correspond to those modes which exited the horizon early on and were, therefore, 
more sensitive to the radiation present and the compact nature of the spatial S s slices. At 
large values of k, the radiation as well as the spatial curvature has inflated away and the 
corresponding modes would fail to notice any sizeable effect due to radiation or curvature. 

With the exact power spectrum in hand, the next question is naturally whether the 
deviations from pure de Sitter are observable. Obviously in short inflation scenarios there 
is a better chance of detecting such an effect. Recall that in closed universes there is a 
minimal wavenumber corresponding to the size of the spatial S* 3 . The longer inflation lasts, 
the higher the wavenumbers are that are responsible for the observed CMB, and thus the 
smaller the effect. Inflation must have lasted long enough to solve the horizon and flatness 
problems, however, and observability thus hinges on the minimal length allowed by these 
constraints. Flatness in particular constrains the present scale factor a p in terms of the 
observed relative density 

The largest scale visible in the CMB, k = a p H p , must therefore be equal to or larger than 

Using the observed value 0, o b s — 1.02, the maximal correction to the power spectrum is 
therefore of order 10~ 2 . A red tilted correction of this size is in practical terms barely 
distinguishable due to cosmic variance. 



5. Discussion and Conclusion 

To summarize, a radiation deformed inflationary phase is an interesting extension to the 
current minimal paradigm. A small amount of radiation is in fact quite generic. Surpris- 
ingly, the radiation deformed geometry allows for an exact analysis of probe scalar field 
perturbations. Calculating the power spectrum of scalar field fluctuations we found that 
the relative power at low multipoles is enhanced. Both the curvature and the radiation 
component contribute at the same order and together give rise to a red power spectrum. 
The observed near critical density of the universe, however, implies that the magnitude 
of this effect is barely detectable. Finally, in the context of a Hartle-Hawking "tunneling- 
from-nothing" scenario, the small amount of radiation expected to be present [9, 10, 26] 
implies that the Euclideanized geometry - a barrel rather than a sphere - does no longer 
single out a unique vacuum state. For the barrel there is no analyticity constraint that 
selects the the Bunch-Davies initial state. 

We conclude with a (partial) list containing some suggestions on how our results might 
connect to other work. 
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• As pointed out, in (short) closed inflation scenarios, the radiation effects are most 
prevalent. In the context of the string theory landscape of de Sitter vacua, a Hartle- 
Hawking "tunneling from nothing" origin of the universe might favor short inflation 
scenarios over inflationary phases with a large number of e-folds because it is far easier 
to arrange a few e-folds in the landscape. Assuming a random form of the de Sitter 
potentials in the landscape, a statistical analysis of the possibility of inflation can be 
done as in [27]. In fact, the proposal in [28] ("living dangerously") suggests that the 
physical parameters of our universe should be on the verge of violating important 
experimental bound (favoring, for example, short inflation). 

• In [3] an initial state of the universe is constructed from closed string tachyons. 
This initial state corresponds to a thermal excitation above the Bunch-Davies initial 
state. Clearly, due to backreaction of the thermally excited state (i.e. radiation), the 
geometry will be modified and one might expect the power spectrum to be affected 
in a similar way as discussed here. 

• In brane inflation scenarios, one expects a period of collision and annihilation of 
branes prior to the last phase of brane inflation directly preceding the reheating 
of our universe. In such scenarios, again assuming that short inflation is favored 
over long periods of inflation, we expect that there will generically be a radiation 
component (due to brane annihilation) during inflation. 

• It is interesting to compare our result with the warm inflation scenario [2]. Warm 
inflation consists of a self interacting inflaton field with a dissipative mechanism 
that generates a sizeable radiation density during inflation. So instead of having a 
reheating stage following a supercooled inflationary phase, radiation is constantly 
generated during inflation. As discussed in [2], for certain models of the inflaton 
potential and dissipation mechanisms there can be a blue, instead of a red, tilt in the 
density perturbation for small wavenumbers due to thermal effects. 

• The authors of [29] consider a scenario where the initial state of the universe is de- 
scribed by a positive cosmological constant and radiation (motivated by the modified 
wave function discussed in [9,10,26]). The subsequent evolution of the universe is 
quite interesting. The presence of a periodic potential (motivated by the string land- 
scape) leads to a Bloch band for the wave function, and the universe cascades down 
this band structure (having originated somewhere near the top of the band as implied 
by the modified wave function) as it evolves. In the process, the energy in the positive 
cosmological constant gets converted into radiation until there is too little energy in 
the cosmological constant to support an inflationary phase. The initial state of a 
scalar field in this scenario will again be described by the mode functions that we 
have found in this paper. The subsequent evolution of the mode functions through 
the cascades warrants a more detailed study. 
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A. Appendix: Heun's and Lame's Differential Equation 

Heun's differential equation is of Fuchsian type with regular singularities at z = 0, 1, a, oo. 
The canonical form of Heun's differential equation is 

d 2 y (1 5 ; e \ dy a(3z 



dz 2 \z z — 1 z — a) dz z(z — l)(z — a)^ ^ ^ 

where y and z are regarded as complex variables and a, j3, 7, 5, e, q, a are complex and 
arbitrary parameters , and a 7^ 0, 1, 00. The first five parameters are linked by the equation 

1 + S + e = a + (3+l (A.2) 

The exponents at the singular points z = 0,1, a, 00 are respectively (0,1 — 7); (0,1 — 
5); (0, 1 — e); (a, (3). According to the general theory of Fuchsian equations the sum of these 
exponents must be equal to 2. This requirement leads to Eq.( |A.2j ), It can be shown that 
any Fuchsian second-order differential equation with four singularities can be reduced to 
the form given by Eq.( [A.l[) which may, therefore, be regarded as the most general form. 

Heun's equation was originally constructed as a generalization of the hyper geometric 
equation. There are three ways in which Heun's equation degenerates to hypergeometric 
equation. 

• Setting a = 1 and q = aft yields the canonical hypergeometric equation for F(a, ^.7; z). 

• Setting e = 0, q = aa(5 also yields the canonical hypergeometric equation for F{a, (3, 7; z). 

• Finally setting a = q = in Eq. (|A.l|) yields the hypergeometric equation for F(a, f3, a+ 
P-6 + l;z). 

Lame's equation is a special case of Heun's general equation (Eq. flA.lj )) for which 

7 = 6 = e = \ , (A.3) 
and thus a + = | . It is common to write 



a = --v, /9 = -(i/ + l), q = --ah. (A.4) 
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With these redefinitions the equation takes the canonical form of Lame's equation 

2/ (*) + n - + 7 + y (*) + J7 7T7 — —Mz) = A.5 

2 \z z — 1 z — aj 4z(z — l)(z — a) 

The parameter v is called the order of the equation. 

Solutions can be constructed to Heun's or Lame's differential equation both as a Frobe- 
nius series in z and as a series in hypergeometric functions. In this paper we have made use 
of the series solution. Upon constructing a series solution to Heun's equation one typically 
obtains three-term recursion relations. It is, in general, impossible to write down the gen- 
eral solution for three-term recursion relations, unlike two-term recursion relations. One 
has to construct the coefficients of the different powers order by order. For further details 
and discussions on Heun's differential equation see reference [21]. 

B. Appendix: The Horizon 

When we work with global coordinates in a de Sitter spacetime, at late times we simply 
get the usual flat slice results. So the horizon is Let us find the late time horizon for 
the geometry at hand. The event horizon at time t is given by 

r°° dt' 

Using Eq.(|2~2|) for a(t), we get 

1(f) = -^==^Ai + co S h\Lt) (lC-,5) - 1(6,5)) (B.2) 



where 6 = arcsin(tanh(Lt)), 5 = A/vI + A 2 , and I(a,(5) is the elliptic integral of the first 
kind defined as 



Jo vl-/3 2 sm^ 



At late times 6 — > ir/2. We shall use the following identities 

7T „ 7T„/1 1 , r2 



1(9,5) = eF {\\^ 5 ') -sin^ cos (a + sin 2 # + sin 4 6 + ..A (B.4) 

where F (h, |; 1; <5 2 ) is the hypergeometric series, and the coefficients are given by 

°» = ir (?? 1;42 ) -1 ' 
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At late times, y A 2 + cosh 2 (Lt) -> e Lt /2, and ir/2-6 -> 2e~ Lt . So we get the following 
late time behavior for the horizon 



where we use the identity F(l/2, 1; 1; S 2 ) = - S 2 and the fact that 6 2 = A 2 /(l + A 2 ). 

So at late times the horizon is the same as the de Sitter horizon L -1 . 
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